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We calculate the admittance of a two-dimensional quantum 
point contact (QPC) using a Boltzman-like kinetic equation 
derived for a partial Wigner distribution function in an ef- 
fective potential. We show that this approach leads to the 
known stepwise behavior of the admittance as a function of 
the gate voltage. The emittance contains both a quantum in- 
ductance determined by the harmonic mean of the velocities 
for the propagating electron modes and a quantum capaci- 
tance determined by the reflected modes. 



PACS numbers: 05.60; 72.10. Bgn; 72.30 

Recent technological progress in manufacturing nano- 
scale solid state structures has made it possible to 
fabricate devices in two-dimensional electronic systems 
(2DES) which operate in the quantum ballistic regime. 
A particular system that has attracted much attention is 
the ballistic quantum point contact (QPC) Jlj-^12]]. 

In this letter we investigate a. c. -electron transport 
through a QPC. Despite detailed investigations, both 
experimental [|]-[| and theoretical pjl|-|l2], of the d.c. 
transport, a.c. transport effects have so far not drawn 
comparable attention, but see recent important papers 
[p|-|TT[ (and |l3|-|l6j for the kinetic response of the reso- 
nant tunnel junction). 

Several experiments |^-|^] have exhibited such quan- 
tum coherent phenomenon as quantization of the d.c- 
conductance versus the gate voltage (or the number of 
propagating modes through the QPC). The theory of 
this phenomenon |iH^-|T2]| explains the d.c. -conductance 
quantization as a consequence of the adiabaticity of the 
electrons' motion through the QPC with smooth bound- 
aries. In an adiabatic geometry (see Fig.l), which is 
smooth on the scale of the Fermi wavelength, the lon- 
gitudinal and transverse motion of electrons can be sep- 
arated in the Schrodinger equation In this case, the 
number of the transverse quantization modes is an adia- 
batic invariant, and the transverse energy plays the role 
of the potential energy for the one-dimensional longitudi- 
nal motion of each mode. Depending on whether the to- 
tal energy of a given electron state is greater or less than 
the effective potential energy of each mode, the mode is 
propagating or non- propagating (see Fig. 2). 

Frequency dependent transport phenomena in QPCs 



are clearly of interest, for the a.c. frequency introduces 
a time-scale that may reveal qualitatively new effects. 

The a.c. -transport has been considered by M. Biittiker 
et al. p-pH, who showed that it is described at low fre- 
quency u>, by an a.c- admittance, Y = l/Z = G — iuj£ 
(Z is the impedance). G is the conductance and £ is 
the emittance, a term first introduced in this context by 
M. Biittiker @. Christen and Biittiker Q derived the 
general expression for the electrochemical capacitance 
and for the displacement current, and established the 
steplike behavior of the QPC emittance in synchronism 
with the conductance steps, and discussed the admit- 
tance for quantized Hall conductors. In paper p-pJ|, the 
omittance was expressed in terms of the geometric ca- 
pacitance, transmission probability, and the densities of 
states of the "mesoscopic capacitor plates'. 

Below we develop a simple method for calculating the 
transport characteristics for the QPC, which is based 
on the Wigner distribution function (WDF) formalism. 
(The effectiveness of the WDF approach to modeling of 
the mesoscopic devices was demonstrated in p^| .)Our ap- 
proach allows us to represent the emittance in terms of 
the capacitance and the inductance, which are expressed 
in the explicit form through the microscopic character- 
istics. Assuming adibaticity, we derive a Boltzmann-likc 
quantum kinetic equation for the partial WDF for the 
QPC. This equation allows us to treat the 2DES in a 
QPC in terms of classical trajectories for the effective ID 
motion. We calculate the a.c. -admittance of the QPC, 
considering propagating and non-propagating (reflected) 
electron modes and show that the real part of the admit- 
tance (the conductance) is quantized versus of the gate 
voltage, consistent with earlier calculations Q based on 
the Landauer formula . We also show that the emit- 
tance £ has a negative part (the quantum inductance), 
which arises from the propagating electron modes and 
whose value is determined by the harmonic mean of the 
electron velocities. The non-propagating electron modes 
give a positive contribution (the quantum capacitance) to 
£ . This contribution is determined by the which depends 
on a geometrical form of the QPC, which is controlled by 
the gate voltage. 

To find the conductivity of a 2DES of a QPC (see 
Fig. 1) taking account of both frequency dependence 
and spatial dispersion, we shall use an approach based 
on the Wigner distribution function |]||o|, /F(f) = 



1 



/ dfTr{pexp[-j-(p+^A(r^)f}9 + {r-r^/2)^(f+f , /2)}. 
Here p is the density matrix operator of the system and 
obeys the quantum Liouville equation; \£ + (r) (9(f)) is 
the Fermi operator creating (annihilating) a particle at 
r; and A is the vector-potential of the electromagnetic 
field. When the characteristic scale of the spatial inho- 
mogeneities exceeds both the radius of interaction among 
the particles and the electron de Broglie wavelength, the 
kinetic equation for the WDF takes the form equivalent 
to the classical kinetic equation |2Cj ]. From the WDF, we 
can find the charge density p and the current density j 

In the adiabatic model of the QPC shown 

in Fig. 1, one assumes the constriction to be 
smoothly varying d'(x) ~ d(x)/L <C 1, where 
2L is a characteristic adiabatic length scale of the 
QPC. With this assumption, the variables in the 
Schrodinger equation can be separated, and the eigen- 
wave function can be written in the form: ip n (x,y) = 
ip n (x)&n[y, d(x)], where the transverse wave function 
*„(!/) = (l/^/d(x))sm{nn[y + d(x)}/2d(x)}9[d 2 (x) - 
y 2 ], satisfies the boundary conditions &n(y)\ y =±d(x) = 0, 
and 9(x) is the Heaviside-step function. The effective 
Hamiltonian for the longitudinal wave function ip n (x) is 
H = — (h 2 /2m)d 2 x + e n (x) + ecj)(x), where cj>(x) is the 
electric potential averaged with respect to the transverse 
coordinate y. 

Adibaticity and the quantization of the transverse mo- 
tion imply that the energy e n (x) in the Hamiltonian H 
has the form e n (x) = -n 2 n 2 Ti 2 /8md 2 (x), with the trans- 
verse quantum number n an adiabatic integral of motion. 
Thus the motion of electrons in the QPC can be viewed 
as a set of n effective one-dimensional electron systems. 
Each effective electron system is influenced by both the 
potential e n (x) and the self-consistent electric potential 
4>(x). To describe these separate systems we introduce a 
partial WDF (PWDF) via 



/„ (x,P* 



J dx' exp (— ip x x'/h) Trp 



(1) 



9+(x-x'/2)9 n (x + x'/2). 
Then we can represent the WDF in the form 

oc °° 
fff(f) = X! f™( x iPx) J dy' exp(-ip y y'/h) 



(2) 



(y-y'/2)<f> n>x (y + y'/2). 

The equation for the PWDF (1) can be derived using the 
Wigner transformation [^9| , po[ : 



Of, 
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dt 



dx 



de n (x) 
dx 



zE(x) 



ir 



dp 



(3) 



where p = p x and E(x) — —d(p(x) /dx. The equation 
(3) is derived in the semi-classical approximation, un- 
der the condition kpL ^> I. In terms of the PWDF 
the nonequilibrium charge density and current den- 
sity can be defined as: p(x,y) = Y^=i Pn(x)$i(y), 
j(x,y)=Y^L 1 jn(x)^l(y), where p n (x) and j n (x) are 
the partial charge and current densities: p n (x) — 

(e/nh) J dp f™(x,p) - f^ {0 \x,p) , and j n (x) = 



(e/irhm) J dppf™(x,p) 



Here, f^^°\x,p) is the equi- 



librium PWDF. The kinetic equation (3) for the PWDF 
has the form as the classical kinetic equation in the effec- 
tive potential e n (x). Thus we can solve this equation by 
the method of characteristics, i.e., by using the classical 
trajectories. 

The collision integral in (3) includes quantum tran- 
sitions [po|| and intermixing of the different electron 
modes (the different PWDF). In this letter, we will 
approximate the collision integral by a single mo- 
mentum relaxation frequency, so that /„ = 



f£(x,p)-f; 



W(0) 



where is the equilib- 

rium PWDF. The equilibrium distribution function 
within the adiabatic approximation is given by: 
f^^(x,p) = n F {[p 2 /2m + e n (x) - p]/T}, where n F is 
the Fermi function with the effective chemical potential 
/i — E n (x) , where p, is the equilibrium chemical potential 
of the 2DEG. The characteristics of kinetic equation (3) 
are the phase trajectories of a one-dimensional motion in 
the potential e n (x) , which is determined from the integral 



of motion, 



the total energy: e — p /2m + e n (x) 



const. We consider a symmetrical QPC, d(x) = d(—x). 
In this case the phase portrait is shown in Fig. 2. The 
separatrix lines (the heavy lines) which pass through the 
hyperbolic point p — 0, x — 0, divide the phase space into 
four regions, within which four sets of phase trajectories 
exist. 

The regions of propagating trajectories (e > e n (0)) oc- 
cupy the regions (see Fig. 2): 1) e > e„(0), p > 0; 2) 
£ > e„(0), p < 0; the regions of non-propagating (reflect- 
ing) trajectories are (e < e„(0)): 3) e < e n (0), x > 0; 
and 4) e < e„(0), x < 0. Within each region, one can find 
the solution of the kinetic equation for the PWDF and 
derive the general formula for the partial charge p n and 
the current densities j n . Here we consider the most in- 
teresting case, when the temperature is very low (T — > 0, 
T <C p) , so that we have a clear separation between open 
channels (e n (0) < p) and closed channels (e„(0) > p). 

Transport through a QPC in principle involves non- 
local (integral) operators in that the charge and cur- 
rent densities at a given point x can influenced by the 
electric field within the whole conductor. Nonetheless, 
it is well known that the static conductance is (to a 
good approximation) specified by the potential difference 
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(bias voltage) between the right and left reservoirs and 
that the detailed electric potential profile does not in- 
fluence it significantly m. This result was derived using 
the Landauer formula |12],|l8| in which the conductance 
is determined by the matrix of the transmission coeffi- 
cients of the electron waves corresponding to different 
propagating modes. We can readily show that this re- 
sult follows from our PWDF approach. In the ballistic 
regime, when L < I, (2L is the distance between the 
reservoirs, I is the mean free path) and at to, v — ► 0, 
we obtain for the open channels j n = (2e 2 /h)V, with 

L 

V = J dxE(x), and for the closed channels j n —0- Hence 

— L 

for the static conductance we get the familiar result |Q: 
G = I/V = 2e 2 Af/h, where TV is the number of the 
open channels: Af—[{2kpd(0)}/iv]; hkp=\/2~rnji. Here 
the brackets [• • •] stand for an integer part. Clearly, the 
static conductance does not depend on the specific form 
of the smooth function, d{x). 

More generally, we can use the formalism of the PWDF 
to calculate the admittance at a finite frequency lj. The 
partial current j n is a function of the longitudinal coor- 
dinate x at u> ^ 0. 

The continuity equation, divj + dp/dt = 0, in the QPC 
at ui allows one to show that the total current, 

X 

hot = J2^=i{jn - iw J dx'p n (x')}, . which includes 

— L 

the electron current 'Y^=\3n(x) and displacement cur- 

X 

rent — iuj X^^Li / dx'p n (x'), is independent of the longi- 

tudinal coordinate x. Within the left reservoir the dis- 
placement current vanishes, so the total current is Itot — 
52?Liin( — L), and the admittance can be determined 
as Y = I tot /V = {l/V)J2n=i3n(-L). In the general 
case, we should determine the field E(x) self-consistently 
within the QPC from Maxwell's equations and afterwards 
calculate the admittance. For the present, we consider 
the long- wavelength approximation, where u* 3> u)L n . 
Here i>* is the typical velocity for electrons in the n th 
channel, and L n characterizes the length of a region for 
that channel. For the propagating modes, L n is the dis- 
tance between the reservoirs (L n ~ 2L) and = v n (0). 
For reflecting modes, L n is the doubled distance between 
the turning points (s n (x n ) = p) and the nearest reser- 
voir [L n ~ 2{L — x n )). Hence the typical velocity in 
this case is: v* L = 2vpy / x n (L — x n )/L, vp — ^/2pjm, 
where 2L is, as above, the characteristic adiabatic length 
scale of the QPC. To calculate the current of the prop- 
agating modes, we can approximate the velocity v n as 
v n (x) ~ v n (0) = u*, while for the reflecting modes: 
v n (x) ~ u*y (|a;| — x n )/(L — x n ). Consistent with the 
symmetry of d(x) we assume that E(x) — E(—x) within 

the QPC and take d(x) = do exp (x/L) 2 . The con- 
tribution of the propagating modes to the total current 



is determined by the bias voltage V and is independent 
of the detailed profile of the electric potential inside the 
QPC, while that of the reflecting modes depends on d(x) 
and E{x). We can write the admittance in the form 
Y = G — iu>£, where as shown above the static conduc- 
tance is given by G = (2e 2 /h)Af and (after some calcula- 
tion) one can show that the emittance £ can be expressed 
as 



£ = -G- 



L 



1R 2 A/+Af 



Here v {o) is the harmonic mean of the velocities v* L in 
the propagating modes: l/v^ = {jj-)J2n=i^/ v n- The 
integer J\f determines the number of reflecting modes: 



TV = 



—TV, 2L being the distance 



(2k F d /n) exp[(i/L) 2 ] 

between the reservoirs. The dimensionless parameter 
value £n reflects the form of the electric field in the region 
(x n , L) filled with the electrons of the n th reflecting chan- 



nel: = (3/2) / dx'(E(x/)/V)y/(x' - x n )/{L - x n ). 

X n 

From their respective expressions one sees immediately 
that the contribution to £ of the reflecting modes is pos- 
itive while that of the propagating modes is negative. In 
this situation we can interpret our results in terms of an 
equivalent circuit in which a capacitance and resistance 
are in parallel followed by an inductance in series. The 
admittance of this circuit is 



Y = G — iuj(C - 



AG 2 



wC<G, wA « -. (5) 
G 



The effective inductance in (5) is given by A = 
(?L/Gv^°\ and the effective capacitance C is given in 
(4). Note, that the formula (5) coincides with the general 
expression for the emittance derived in [fll| (see formula 
(7) in [jlH), where the emittance was expressed in terms 
of the geometric capacitance, transmission probability, 
and the densities of states of the "mesoscopic capacitor 
plates" . Our approach allowed us to represent the emit- 
tance in terms of the capacitance and the inductance, 
which are represented in the explicit form through the 
microscopic characteristics. 

The dependence of the emittance on the number of 
open channals shows immediately that it is a stepwise 
function of the gate voltage. When the gate voltage ap- 
proaches a point for which 2kpdj r K is an integer, so that 
an additional channel opens (or closes), our expressions 
for the inductance and the capacitance diverge. In this 
case, our assumption that v* 3> toL n is violated, and the 
contribution of these points into the admittance should 
be calculated separately, using the full self-consistent 
Maxwell equation approach j n = (/x — e ?l (0))/e n (0), 

In conclusion, we have developed a new approach, 
based on a partial Wigner distribution function, to ana- 
lyze the a.c. electron transport properties of a quantum 
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point contact. Treating the ballistic QPC in the adia- 
batic approximation, we derived a Boltzman-like equa- 
tion for the partial Wigner distribution function in an 
effective potential brought about by the quantized trans- 
verse modes. We analyzed this equation in terms of prop- 
agating and reflecting trajectories in the semi-classical 
approximation. 

Our results establish that the a.c. electron transport 
depends directly on the the number of propagating and 
reflecting modes and that certain features are sensitive 
to the form of the distribution of the electric field in 
the QPC. In particular, the real part of the admittance 
(the conductance) is determined by the number of prop- 
agating electron modes and does not depend on the spa- 
tial distribution of the electric field inside the QPC |Q. 
The imaginary part of the admittance (the emittance) 
exhibits stepwise behavior as a function of the gate volt- 
age and consists of two parts: the quantum inductance 
and the quantum capacitance. The quantum inductance 
is determined by the harmonic mean of the velocities for 
the propagating electron modes. The quantum capaci- 
tance is specified by the reflecting modes and is sensitive 
to the geometry of the QPC. The emittance can be con- 
trolled by the gate voltage 

We stress that the effective quantum inductance and 
capacitance, and the equivalent circuit, are concepts valid 
within the linear response, low-frequency approximation. 

For the high-frequency case, and when new propagat- 
ing and non-propagating modes can appear or disappear, 
the frequency dispersion of the admittance is more com- 
plicated than the linear one given by the equivalent cir- 
cuit of Eq. (5) . This case must be considered using the 
self-consistent Maxwell equations for the electric field in 
the QPC. We are presently investigating this problem. 
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FIGURE CAPTIONS 



Fig. 1: The geometry of the microconstriction. The 
width is denoted by 2d{x), the narrowest width is 2do, 
and the effective length is 2L. 

Fig. 2: The plane of phase trajectories for one- 
dimensional motion determined by conservation of inte- 
grals of motion. The heavy lines are separatrices that 
separate the propagating modes (sections 1 and 2) and 
non-propagating (reflecting) modes (sections 3 and 4). 
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